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Abstract. In this note, we summarize and compare various model predictions for pp 

r\] • total cross-section a^^^ , giving an estimate of the range of predictions for the total cross- 

^^ ' section, af^^ expected at the LHC. We concentrate on the results for a\^^ obtained in a 

f^ I particular QCD based model of the energy dependence of the total cross-section, including 

v^ , the effect of soft gluon radiation. We obtain the range of predictions in this model by 

is^ , exploring the allowed range of model parameters. We further give a handy parametrisation 

^^ . of these results which incidentally spans the range of various other available predictions 

fH , at the LHC as well. 
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1. Introduction 



Oh! 
> 

X 
^. 

5t 1 Energy dependence of total hadronic cross-sections has been the focus of intense 

theoretical interest as a sensitive probe of strong interactions long before the es- 
tablishment of QCD as 'the' theory of hadrons. Even now, notwithstanding cred- 
itable successes of perturbative and lattice QCD, a first principle description of 
total/elastic and inelastic hadronic cross-section is unavailable. More pragmati- 
cally, for a correct projection of the expected underlying activity at LHC, a reliable 
prediction of total non-diffractive cross-section is essential to ensure the extraction 
of new physics from the LHC data. Surely we will have to depend - at the initial 
stages of LHC - upon predictions based on our current understanding of these mat- 
ters. Only much later it may become feasible to use the LHC data itself towards 
this goal. Hence, a critical evaluation of the range of theoretical predictions, is 
absolutely essential. 

The hadronic cross-section data exhibit, and require explanation of, three basic 
features: 
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(i) the normalization of the cross section, 

(ii) an initial decrease and 

(iii) a subsequent rise with energy. 

Various theoretical models exist which are motivated by our theoretical understand- 
ing of the strong interactions. The parameters in these models, in most cases, are 
fitted to explain the observed low energy data and the model predictions are then 
extrapolated to give the a^^^ at the LHC energies. There are different classes of 
models. The highly successful Donnachie-Landshoff (DL) parametrisation [1] of 
the form 

atot(s)=Xs' + Ys-\ (1) 

has been used for a very long time. Here the two terms are understood as arising 
from the Regge and the Pomeron trajectories, the e being approximately close to 
zero and 77 being close to 0.5. These values seem to be consistent with a large, 
but not all, body of the hadronic cross-sections. In this note we will first present 
phenomenological arguments for the approximate values of these parameters which 
seem to be required to describe the data satisfactorily. As a matter of fact, there 
also exist in the literature discussions of the 'hard' Pomeron [2] motivated by the 
discrepancies in the rate of energy rise observed by E710 [3], E811 [4] and the 
CDF [5]. In addition, a variety of models exist wherein the observed energy depen- 
dence of the cross-section, along with few very general requirements of factorisation, 
unitarity and/or ideas of finite energy sum rules (FESR) , is used to determine the 
values of model parameters [6-12]. The so-obtained parametrisations are then ex- 
tended to make predictions at the LHC energies. There also exist QCD motivated 
models based on the mini-jet formalism [13-15], wherein the energy rise of the total 
cross-sections is driven by the increasing number of the low- a; gluon-gluon collisions. 
These models need to be embedded in an eikonal formalism [16] to soften the violent 
energy rise of the mini-jet cross-sections. Even after eikonalisation the predicted 
energy rise is harder than the gentle one observed experimentally [15,17]. A QCD 
based model where the rise is further tamed by the phenomenon of increasing emis- 
sion of soft gluons by the valence quarks in the colliding hadrons, with increasing 
energy [18,19], offers a consistent description of crfot. Thus we have a variety of 
model predictions for crf^j at the LHC. In this note we compare these predictions 
with each other to obtain an estimate of the 'theoretical' uncertainty in them. 



2. Phenomenological models 

The two terms of eq. (1) [1] reflect the well-known duality between resonance and 
Regge pole exchange on the one hand and background and Pomeron exchange on 
the other, established in the late 60's through FESR [20]. This correspondence 
meant that, while at low energy the cross-section could be written as due to a 
background term and a sum of resonances, at higher energy it could be written as 
a sum of Regge trajectory exchanges and a Pomeron exchange. 

It is well to ask (i) where the 'two-component' structure of eq. (1) comes from and 
(ii) why the difference in the two powers (in s) is approximately a half. Our present 
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knowledge of QCD and its employment for a description of hadronic phenomena 
can be used to provide some insight into the 'two-component' structure of eq. (1). 
This begins with considerations about the bound state nature of hadrons which 
necessarily transcends perturbative QCD. For hadrons made of light quarks (g) 
and gluons (g), the two terms arise from qq and gg excitations. For these, the 
energy is given by a sum of three terms: (i) the rotational energy, (ii) the Coulomb 
energy and (iii) the 'confining' energy. If we accept the Wilson area conjecture in 
QCD, (iii) reduces to the linear potential [21,22]. Then the hadronic rest mass for 
a state of angular momentum J can be obtained by minimising the expression for 
the energy of two massless particles (qq or gg) separated by a distance r. 

Explicitly, in the CM frame of two massless particles, either a gg or a gg pair 
separated by a relative distance r with relative angular momentum J, the energy 
is given by 

E,{J,r) = — -^ + arr, (2) 

r r 

where i = 1 refers to qq, i = 2 refers to gg, r is the 'string tension' and the Casimir's 
are Ci — Cp = 4/3, C2 — Cq = 3. a is the QCD coupling constant evaluated at 
some average value of r and whose precise value will disappear in the ratio to be 
considered. The hadronic rest mass for a state of angular momentum J is then 
computed through minimising the above energy 



Mj(J) =min 

r 

This is then given by 



2J C\a 
r r 



Cr 



(3) 



M,{J) = 2^iaT)[2J-aa]. (4) 

This can then be used to obtain the two sets of linear Regge trajectories 

"^(■5) == —k~ + ( Z7^~ ] * = Q;j(0) + a'^s, (5) 



2 \8C, 

Note that ai are not the coupling constants. 
Thus, the ratio of the intercepts is given by 



«gg(0) 
"^(0) 



Cg/Cf = 7. (6) 



Employing our present understanding that resonances are qq bound states while 
the background, dual to the Pomeron, is provided by gluon-gluon exchanges [22], 
the above equation can be rewritten as 

^ - Co/C. = l- (7) 

o:r(Oj 4 

If we restrict our attention to the leading Regge trajectory, namely the degenerate 
p — uj — (j) trajectory, then aR(0) = 77 w 0.48-0.5, and we obtain for e « 0.08-0.12, 
a rather satisfactory value. The same argument for the slopes gives 
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Hence, if we take for the Regge sfope a^ « 0.88-0.90, we get for ap « 0.39-0.40, 
in fair agreement with lattice estimates [23]. 

We now have good reasons for a break up of the amphtude into two components. 
To proceed further, it is necessary to reahze that precisely because massless hadrons 
do not exist, eq. (1) violates the Froissart bound and thus must be unitarized. To 
begin this task, let us first rewrite cq. (1) by putting in the 'correct' dimensions 

atot(s)='Ti(s/s)^ + a2(s/s)i/2, (9) 

where we have imposed the nominal value r\ = 1/2. It is possible to obtain [19] 
rough estimates for the size of the parameters in eq. (9). A minimum occurs in 
^tot(s) at s = s, for 02 = 2ea'i. If we make this choice, then eq. (9) has one less 
parameter and it reduces to 

atot(s)=ai[(s/sr + 2e(s/s)i/2]. (10) 

We can isolate the rising part of the cross-section by rewriting the above as 

Stot(s) - ai[l + 2e(s/s)i/2] + ^i[(s/s)^ - 1]. (11) 

Equation (11) separates cleanly the cross-section into two parts: the first part is 
a 'soft' piece which shows a saturation to a constant value (but which contains no 
rise) and the second a 'hard' piece which has all the rise. Moreover, s naturally 
provides the scale beyond which the cross-sections would begin to rise. Thus, our 
'Born' term assumes the generic form 



-^tot 



(s) = Crsoft(s) + "d{s - s)crhard(s) (12) 



with CTsoft containing a constant (the 'old' Pomcron with ap(0) ~ 1) plus a 
(Regge) term decreasing as l/\/s and with an estimate for their relative magni- 
tudes (a^jox ~ 2e). In the eikonalised mini-jet model used by us [19] the rising 
part of the cross-section dhard is provided by jets which are calculable by pertur- 
bative QCD, obviating (at least in principle) the need of an arbitrary parameter e. 
An estimate of a\ can also be obtained [19] and is ^ 40 mb. 

As said earlier, the DL parametrisation [1] is a fit to the existing data of the 
form given by eq. (1), with e = 0.0808, 77 ~ 0.4525. This fit has been extended to 
include a 'hard' Pomeron [2] due to the discrepancy between different data sets. 
The BH model [6] gives a fit to the data using duality constraints. The BH fit for 
cr* = crPP/crPP as a function of beam energy v^ is given as 



o-'*' = Co + ci 



ln(i^/m) + C2ln^(i^/m) +/3p/(i^/to)^ ^ ±8(vlray- 



where /i — 0.5, a = 0.415 and all the other parameters in mb arc cq = 37.32, c\ — 
-1.440 ± 0.07, C2 = 0.2817 ± 0.0064, /3p. = 37.10, ,5 = -28.56. The fit obtained by 
Igi and Ishida [8] using the finite energy sum rules (FESR) gives LHC predictions 
very similar to those given by the BH [6] fit. Avila et al [9] gave a fit using 
analyticity arguments whereas Cudell et al [10] gave predictions at the LHC energies 
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by extrapolating fits obtained to the current data again using constraints from 
unitarity, analyticity of the S-matrix, factorisation, coupled with a requirement 
that the cross-section asymptotically goes to a constant as ln(s) or In^(.s), in the 
framework of the COMPETE program. 

In the mini-jet models the energy rise of a^^^ is driven by the increase with energy 
of the CTj^f given by 



ryS/2 .1 r-1 

^j^f (■«) = / dpt / dxi / dx2 2Z fi\A{xi)f:j\B{x2) 

x%:M(f), (13) 

dpt 

where subscripts A and B denote particles (7,^, . . .), i,j,k,l are partons and xi,a;2 
the fractions of the parent particle momentum carried by the parton. s — X1X2S 
and (T are hard partonic scattering cross-sections. This parton model used in the 
calculation is illustrated in figure 1. Note here that the experimentally measured 
parton densities in the proton and the elements of perturbative QCD is all the 
input needed for the calculation of cr^^ . The rate of rise with energy of this cross- 
section is determined by pt min and the low- a; behaviour of the parton densities. As 
said before, the rise with energy of this cross-section is much steeper than can be 
tolerated with the Froissart bound. Hence it has to be embedded in an eikonal 
formulation given by 

af£ = 2 fd^l - c-3"x-^"(''^-)], (14) 

where Rex'^-^(6, s) = and 23to x"^^(6, s) = n^^{b,s) is the average number 
of multiple collisions which are Poisson distributed. As outlined in eq. (11) this 
quantity too has contributions coming from soft and hard physics and can be written 
as 

n-4^(6,,s)=nf,f+<f, 

^ Aii{b)at£,{s) + Af,f(b)af,^{s). (15) 

In the second step the number n(6, s) has been assumed to be factorisable into an 
overlap function A{b) and the cross-section a. 

The sketch in figure 2 indicates the relationship between the multiple scatterings 
and the overlap function. The assumption of factorisation as well as the split up 
between the two contributions, hard and soft, for the n(6, s) are only approximate. 
The extent to which this softens the energy rise, depends on the b dependence of 
n{b,s), i.e., that of A{b) in the factorised case. The normal assumption of using 
the same form of A{b) for both the hard and the soft part, given by the Fourier 
transform of the electromagnetic form factor (FF), still gives too steep a rise even 
in this eikonalised mini-jet model (EMM) [16]. In our model this rise is tamed 
by including the effect on the transverse momentum distribution of the partons in 
the proton, of the soft gluon emission from the valence quarks in the proton [18]; 
the effect increasing with increasing energy. In this description, the transverse 
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A + B * jet + jet + X 



Figure 1. Parton model picture for jet production. 




Figure 

hadrons. 



2. Multiple scattering giving rise to the transverse overlap of the 



momentum distribution of the quarks in a proton can then be calculated in a semi- 
classical picture as arising from the re-summation of a large number of soft gluons. 
This in turn allows us to calculate the transverse overlap function. 

Figure 3 sketches the multiple emissions of gluons which gives rise to the trans- 
verse momentum distribution of the valence quarks and hence the overlap function. 
The non-perturbative soft part of the eikonal includes only limited low energy gluon 
emission and leads to the initial decrease in the proton-proton cross-section. On 
the other hand, the rapid rise in the hard, perturbative jet part of the eikonal is 
tamed into the experimentally observed mild increase by soft gluon radiation whose 
maximum energy (^max) rises slowly with energy. Thus the overlap functions A{b) 
are no longer a function of b alone. We denote the corresponding overlap function by 
^BAf(&, 9max)- The CT^^j Can then be computed using eq. (14) where 2'^mx^^{b,s) 
is given by 

2^mx{b,s) = n(6, S;ginax,Ptmin) = "soft +njot 

= Asoltib, s)crsoft(s) + AsNib, gmax)crjot(s). (16) 
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valence quark 



Figure 3. Multiple gluon emission giving rise to the transverse momentum 
distribution of the valence quarks. 



The function Abjv(6, (/max) [18] is determined by gmax, which in turn depends on the 
energy and the kinematics of the subprocess. What we use is its average value over 
all the momentum fractions of the parent partons. We need to further make a model 
for the 'soft' part Asott, which is determined by the non-perturbative dynamics. It 
is this part of the eikonal that contributes to the a^^^ at high energies through its 
impact on the turn around from the decreasing Regge behaviour to the softly rising 
behaviour around ^/s ~ 15 GeV, where the hard part contribution is minuscule. 
Note that we have taken Usoit to be factorised into a constant soft cross-section asoh 
and taken Asoh to be given by the function ABN{b,q^^x) except for the fact that 
it is not possible to calculate the (/max ^'^'^ the soft processes, as in the case of hard 
processes. We further postulate that the <7max is the same for the hard and soft 
processes at low energy, parting company around 10 GeV where the hard processes 
start becoming important. A good fit to the data requires (/max at low energies to 
be a very slowly decreasing function of energy, with a value around 0.20 MeV at 
^/s = 5 GeV rising to about 0.24 MeV, -^/s > 10 GeV, the upper value of this soft 
scale being completely consistent with our picture of the proton. Further, we need 
to fix one more parameter for non-perturbative region, the Usoft- For the pp case 
it is a constant ag which will fix the normalization of a^^^ , whereas for the pp the 
duality arguments suggest that there is an additional -y/s dependent piece ~ l/\/s- 
Thus neglecting the real part of the eikonal, n{b, s) in our model is given by 

n{b, S) = ABNib, CfxXoT + ^BN{b, CL)'^jet(s;Ptmi„), (17) 

where 

<!h--o, aZ,-<y.(l + ^)- (18) 

Thus the parameters of the model are pt min and ctq . In addition, the evaluation 
of Abn involves as in the infrared region, for which we use a phenomenological 
form inspired by the Richardson potential [18]. This involves a parameter p which 
for the Richardson potential takes value 1. Values of pt min, era and p which give 
a good fit to the data with the GRV parametrisation of the proton densities [24] 
are 1.15 GeV, 48 mb and 3/4 respectively, as presented in ref. [19]. These values 
are consistent with the expectations of the general argument [19]. We expect these 
best fit values to change somewhat with the choice of parton density functions 
(PDF). Since we are ultimately interested in the predictions of the model at TeV 
energies, we need PDF parametrisation which cover both the small and large Q^ 
range (2 < Q^ < 10**) as well as are valid up to rather small values of x(~ 10~^). 
Further, since our calculation here is only LO, for consistency we have to use LO 
densities. We have repeated the exercise then for a range of PDF's [25-27] meeting 
these requirements. For each PDF, it is onset of the rise that fixes the pt min, 
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GRV940LD densities 

Blocl<-Hal2en model 

-- a„=48mbp=0.75p„„=1.15GeV 




proton-proton 

A proton-antiproton 

• UA5 
A UA1 
o UA4 

* CDF 
D E710 
T ES11 



10 10 

V's ( GeV ) 




GRV98LO densities 
— Block-Halzen model 

- ci„=50mbp=0.77p,„i„=1.10 
-- <j„=45mbp=0.70p,„i„=1.10 




proton-proton 

A proton-antiproton 

• UA5 
A UA1 
o UA4 

* CDF 
n E710 
' E811 



10 

•is { GeV ) 



Figure 4. Comparison of tlie G.G.P.S. model [19] predictions for GRV [24] 
and GRV981o [26] densities witli tlie BH [6] predictions. 



ctq controlling the normalisation and p determining the slope of the rising part of 
the cross-section. We find that it is possible to get a satisfactory description of 
all the current data, for all the choices of PDF's considered. The corresponding 
range of values of ptmin, ctq and p are given in table 1. Figure 4 compares with the 
data the predictions of our QCD based model along with those obtained by the 
considerations of unitarity and factorisation [6]. As can be seen clearly, both are 
able to describe the current data on total cross-sections equally satisfactorily. As 
can be seen from figures 4 and 5, it is possible to get equally satisfactory description 
of the data in our QCD based model, for all the chosen PDF's, by tuning the soft 
parameters by a small amount. In the next section we compare these model results 
with the predictions of all the other models. 



3. Model predictions for cr^^ at the LHC 

Figure 6 summarises the predictions of the different models described in the previous 
section. The shaded area gives the range of predictions in the eikonalised mini- 
jet model with soft gluon re-summation [19] (the G.G.P.S. model), the different 
PDF's used giving the range as described in the earlier section. The solid line gives 
prediction obtained using the GRV parton densities [24] in the model. The curve 
d indicates predictions of the DL fit [1]. The (BH) curve c and the uppermost 
curve a, are the results of two analytical models incorporating constraints from 
unitarity and analyticity, from [6] and [9], respectively. The predictions obtained 
by Igi and Ishida, using FESR follows very closely to that given by the BH curve. 
Further, curve b is the result of a fit by the COMPETE collaboration [10]. The 
parametrisation for the DL curve and BH curve is already given in the last section. 
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Figure 5. Predictions of the G.G.P.S model [19] for the GRV94Io [25] and 
the MRST [27] densities for model parameters mentioned in table 1. 



It is gratifying to see that the range of results of our QCD motivated mini-jet 
models for the LHC span the other predictions obtained in models using unitarity, 
factorisation, analyticity along with fits to the current data. Thus the two sets of 
predictions seem consistent with each other. 

We have parametrised the results of our EMM model with a In (s) fit. We found 
that in most cases this gave a better representation of our results than a fit of the 
Regge-Pomeron type of the form of eq. (1). The top edge of the EMM model 
prediction is obtained for the MRST parametrisation whereas the lower edge for 
the GRV981o. We give fits to our results for a^^ of the form 



Table 1. Values of ptmin and uo corresponding to the dif- 
ferent parton densities in the proton, for which the EMM 
(as described in ref. [19]) gives a satisfactory description of 



PDF 


Ptmin (GeV) 


(To (mb) 


P 


GRV [24] 


1.15 


48 


0.75 


GRV9410 [25] 


1.10 


46 


0.72 




1.10 


51 


0.78 


GRV9810 [26] 


1.10 


45 


0.70 




1.10 


50 


0.77 


MRST [27] 


1.25 


47.5 


0.74 




1.25 


44 


0.66 
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G.G.P.S. model, Phys. Rev. D 72, 076001 (2005) 

— G.G.P.S. model, using GRV P.D.F. 

modified G.G.P.S. model, using GRV P.D.F. 
Luna-Menon model, hep-ph/oi 05076 
Gudell et. al. model, hep-ph/0212101 
Block-Halzen model, Phys.Rev. d 72 036006 (2005) 
Donnachie-Landshoff model, prl b296(1992) 227,Yy 




proton-proton 

proton-antiproton 

UA5 

UA1 

UA4 

CDF 

E710 

E811 



10 



10 



10 



10 

Vs ( GeV ) 



Figure 6. Predictions for a^^^ in various models. The shaded area 
gives the range of results in the eikonalised mini-jet model with soft gluon 
re-summation [19] (the G.G.P.S. model) the solid line giving the prediction 
obtained using the GRV parton densities [24] in the model. Curve d indi- 
cates predictions for the DL fit [1]. Curve c and the uppermost curve a, are 
the results of two analytical models incorporating constraints from unitarity 
and analyticity, from [6] and [9], respectively. The prediction obtained by Igi 
and Ishida, using FESR follows very closely to that given by the BH curve. 
Further, curve b is the resuH of a fit by the COMPETE collaboration [10]. 



ao + ais'' + 02 ln(s) + a^ In^(s). 



(19) 



The values of various parameters for the top end lower edge as well as the central 
curve are given in table 2. 
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Table 2. Values of oq, oi, 02, 03 and b parton densities in the proton, for which 

-'tot- 



the EMM (as described in ref. [19]) gives a satisfactory description of a^^ 



ao (mb) ai (mb) b 02 (mb) a-j (mb) 

Top edge 23.61 54.62 -0.52 1.15 0.17 

Center -139.80 193.89 -0.11 13.98 -0.14 

Lower edge -68.73 125.80 -0.16 11.05 -0.16 



4. Conclusions 

We thus see that the range of the results for the o-f^j from our QCD motivated EMM 
model [19] spans the range of predictions made using the current data and general 
arguments of unitarity and/or factorisation. Further, we give In (s) parametrisation 
of the model results for a^^^ which may be used in evaluating the range of predictions 
for the underlying event at the LHC. 
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